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ABSTRACT 

We analyze potential effects of an extraterrestrial civilization’s use of orbiting mirrors to illuminate 
the dark side of a synchronously rotating planet on planetary transit light curves. Previous efforts 
to detect civilizations based on side effects of planetary-scale engineering have focused on structures 
affecting the host star output (e.g. Dyson spheres). However, younger civilizations are likely to be 
less advanced in their engineering efforts, yet still capable of sending small spacecraft into orbit. Since 
M dwarfs are the most common type of star in the solar neighborhood, it seems plausible that many 
of the nearest habitable planets orbit dim, low-mass M stars, and will be in synchronous rotation. 
Logically, a civilization evolving on such a planet may be inspired to illuminate their planet’s dark 
side by placing a single large mirror at the L2 Lagrangian point, or launching a fleet of small thin 
mirrors into planetary orbit. We briefly examine the requirements and engineering challenges of such 
a collection of orbiting mirrors, then explore their impact on transit light curves. We incorporate 
stellar limb darkening and model a simplistic mirror fleet’s effects for transits of Earth-like {R — 0.5 
to 2 i?Earth) planets which would be synchronously rotating for orbits within the habitable zone of 
their host star. Although such an installation is undetectable in Kepler data, JWST will provide the 
sensitivity necessary to detect a fleet of mirrors orbiting Earth-like habitable planets around nearby 
stars. 

Keywords: extraterrestrial intelligence - astrobiology - planets and satellites: terrestrial planets 


1. INTRODUCTION 



search for extraterrestrial intelligence (SETI). Most prior 
searches concentrated on detecting electromagnetic com- 
munications in radio/microwaves and optical/IR bands 
dVon Korff et al.||2013| iKorpela et al.||2011| iHoward et 
| al.|lj007||'i'arier|ljmil |ij' owyer et al.|19i)7[|bagan & Drake 
|l975p. Such searches have operated intermittently since 
the 1960s and continue today with greatly enhanced fre- 
q uency coverage and detection sensitivity. 


Freitas (1983) suggested that the first evidence we find 


of an extraterrestrial civilization may be the detection of 
a physical artifact rather than an electromagnetic signal. 
In this paper we investigate the potential for detection 
of a large fleet of small spacecraft in orbit around a ter¬ 
restrial exoplanet, for example, mirrors illuminating the 
dark side of a planet in synchronous rotation. Such a 
fleet could be observable through its effect on the transit 
light curve. 

The remainder of the introduction motivates this sce¬ 
nario and places it in context with conjectures regard¬ 
ing other extraterrestrial engineering projects and our 
curr ent understanding of exoplanet habitability. In Sec- 
tion |2.1| we explore plausible methods of dark-side illumi¬ 
nation, discussing as socia ted technical and engineering 
challenges in Section |2.2[ Section describes a simple 
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model of the effects of an orbiting mirror fleet on transit 
light curves. We present the resulting light curves in Sec¬ 
tion 1^ focusing on the signature and magnitude of the 
mirrors’ influence, and examining the effects of altering 
various stellar, planetary, and mirror fleet parameters. 
We consider prospects for detection in Section]^ 

1.1. Motivation 

The most common type of star in our Galaxy is the M 
dwarf. These stars have very long lives. Because of their 
low luminosity, the habitable zone (HZ) around M dwarfs 
is narrow and close to the star. Due to tidal considera¬ 
tions, such planets will probably alwa ys present one face 
to the star, except in unusual cases. ( Tarter et al.||2007 l 
Although young ( <100 Myr) dM stars have active chro- 
mospheres and high UV output that could be damaging 
to biologi cal organisms, suc h activity decreases rapidly 


with age (lEngle et ah] 20091) . There is also a wide vari- 

II 


ation in activity at all ages (Shkolnik fc Barman||2014) 


and its precise impact on exoplanet atmospheres, water 


active research (e.g. Mulders et al. 

2014 

Tian||2014 

Zu-1 

luaga et al.|2013 Crentell et al. 12012). bince the specific 


not well-defined, and given that M and K dwarfs are 8 0% 
of stars in the solar neighborhood |Henry et al.||2006 ), it 
seems plausible that most nearby habitable planets are 


plausible that most nearby habitable planets 
in synchronous rotation in orbits around low-mass stars. 
A significant fraction of M dwarfs do host small planets 
within their habitable zones, with the detailed estimates 


depending on the choice of HZ (Dressing & Charbonneau 
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2013 

Gaidos 

2013 |Kopparapu 

2013 

). rent terrestrial civilization 4x10^® erg s ^). Given 

Lite arounc 

an M dwarf coul 

d pose unique challenges the possibility ot extinction, it seems likely that Kar- 


for life, and for the evolution of intelligence. But once 
a civilization has either evolved in planeta or migrated 
from elsewhere, it might wish to utilize portions of the 
planet that do not receive illumination from the star. Af¬ 
ter careful consideration of the climatic effects, the first 
deliberate planetary scale engineering project might be 
illumination of the planet’s dark side. The most straight¬ 
forward method is via orbiting thin mirrors, either plac¬ 
ing a single large mirror at the L2 Lagrange point or, 
more likely, launching a fleet of smaller mirrors into lower 
orbits. To provide significant illumination, these mirrors 
must have a total reflective surface area approximately 
equal to the disk area of the planet. This could render 
them detectable in a planetary transit. 

If the primary motivation is dark-side illumination, 
we anticipate such structures would be found primarily 
around synchronously rotating planets, most likely or¬ 
biting M stars. As we discuss in Section [TT| habitable 
worlds in certain orbits around earlier-type stars (e.g. G 
stars) might also be in synchronous rotation at the cur¬ 
rent time. In addition, an alien civilization might wish to 
utilize a previously uninhabited planet in their star sys¬ 
tem, or simply ensure illumination of the night side of a 
planet in non-synchronous rotation. Other possible uses 
of fleets of comparable area include photoelectric power 
generation. Thus, although in what follows we focus our 
attention on worlds with a perpetual dark side, such en¬ 
gineering projects might occur in other situations, and 
around many types of stars. 

1.2. Prior proposed megastructures 


Dyson (1960) suggested that an extraterrestrial civi¬ 
lization might build a giant structure to capture and uti¬ 
lize most or all of the visible emissions of a star, radiating 
only waste heat. A Dyson sphere would be detectable 
as a stellar-luminosity ^ 300 K featureless blackbody 
source if the structure completely encases the star, or as 
a star with a significant infrared (IR) excess in the case of 
a partial sphere. If civilizations capable of such feats of 
engineering were common, we could likely detect them in 
IR sur veys. Searches of the IR sky sensitive to both com¬ 
plete (Carrigan 2009) and partial (Jugaku & Nishimura 
2004; Conroy & Werthimer 2003, unpublished^; Globus 
et al. 2003, unpublished^ ) Dyson spheres have been 
conducted using the IRAS dataset. While these have 
detected a few interesting candidates, none have been 
shown to indicate intelligent origin. Two high-sensitivity 
searches for Dyson spheres are ongoing, one searching 
for IR signatures of Dyson spheres in the WISE survey 
( Wright et al.|[^14a|b Griffith et al.||2015 ), and another 
examining Kepler Wghi curves tor umisuaf tr ansit eclipses 


that could indicate a partial Dyson sphere (Marcy et al. 

"A civilization advanced enough to trap most of the 
available output of their host star 4x10^^ erg s 


-I'l 


meets the definition of Kardashev Type II on |Kardashev| 
(1964)’s scale of technological advancement of a civiliza- 
tion. A Type I civilization uses energy on the scale of cur- 

^ From http://seti.berkeley.edu/IR_Excess_Search retrieved 09 
July, 2015 


dashev Type II civilizations are much rarer than those 
capable of sending small spacecraft into orbit, and that 
a continuum of possible engineering stages should exist. 
We might consider a civilization with significant ability 
to modify its planetary environment, but not capable of 
b uilding A U. sca le structures as Type 1.5. 


^ Mentioned in Carrigan (2009 I 


Carrigan (2012) recently summarized several other po- 
tential archaeological signatures for civilizations at vari¬ 
ous levels on the Kardashev scale, concluding that apart 
from electromagnetic communications, these are gen- 
erally undetec table with current human technologies. 
Arnold (2005|) suggests that a extraterrestrial civiliza¬ 
tion could deliberately create planetary scale structures 
with unique shapes or periodicities that could be de¬ 
tected when the structure transits the star, but could be 
distinguished from a planetary transit with current tech¬ 
nology. Our proposed engineering project is of similar 
scale, but with a useful function besides communication 
and with a form that follows from its function. 

2. DARK-SIDE ILLUMINATION UTILIZING MIRRORS 
2.1. Methods of dark-side illumination 

There are multiple methods by which an alien civiliza¬ 
tion could illuminate the dark side of a synchronously 
rotating planet. Conceptually, the simplest is placing 
a single large mirror at the L2 Lagrange point (Figure 
[^). To provide illumination equivalent to the bright- 
side stellar illumination, it would need to be a large disk 
or annulus somewhat larger than the planetary diameter. 
Since the mirror is locked in the anti-stellar direction it 
would provide non-uniform directional illumination on 
the ground, with regions near the terminator receiving 
illumination at higher incidence angles than regions near 
the anti-stellar point. One notable difficulty is maintain¬ 
ing positional and directional stability at the unstable L2 
point. A structure and control system that could allow 
pointing and positional control of a 6000 km to 25000 
km diameter mirror would be a significant engineering 
challenge. 

In the case of a single large {Rm > Rp) mirror at the 
L2 Lagrange point, the optical transit light curve will 
last longer and be deeper than if only the planet blocked 
the starlight. If the mirror has been engineered so all 
intercepted light is redirected to the planet, and none 
towards our line-of-sight, the planetary eclipse phase will 
be the same as if the mirror were not present. Thus 
the difference in transit and eclipse time-span would be 
the main clue to the presence of such a structure, but a 
similar timing effect could occur for a planet in a highly 
elliptical orbit with periastron at the time of eclipse and 
apastron at the time of transit. 

A less technically challenging method would be to 
deploy a fleet of independently steerable smaller mir¬ 
rors. Figure schematically illustrates this possibility 
in cross-section for mirrors in a circular orbit. Elongated 
mirror orbits, as depicted in Figure [TJ:, might be prefer¬ 
able, since they spend a greater fraction of their orbital 
period in a position to usefully illuminate the planet’s 
dark side. 

For a fleet of many individual reflectors, the translu¬ 
cent fleet of mirrors would block some starlight, deepen- 
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Figure 1. Schematic illustration of three methods of dark-side 
illumination (not to scale). Planetary gray scale bands indicate 
different levels in stellar illumination. In the three cross-sectional 
drawings, (a) shows a large circular or annular mirror stationed at 
the L2 Lagrange point, (b) shows multiple small mirrors in circular 
orbits, (c) shows multiple small mirrors in elliptical orbits designed 
to maximize the duty cycle of the mirrors. 

ing and lengthening the transit. Before modeling in de¬ 
tail the effects of a fleet of space-based mirrors on transit 
light curves, we first discuss technical and engineering 
considerations involved in its implementation. 


2.2. Technical & Engineering Considerations 

Current materials technologies allow construction of 
metalized plastic films with surface mass densities of 
0.7— 1.5 mg cm“^; neglecting supporting structures, the 
mass of a 1 km^ mirror would be 7000 to 15000 kg 
(Bryant et al.|2014|. F uture materia l s deve lopment (car- 
bon nanotube fabrics: Li ma et al. ( 2011 ), space-based 

WrightI (|iy92| p.l09)) could 
ace mass densities from 0.01 


fabrication of metal films: 
result in materials with sur: 
mg cm“^ to as little as 0.002 mg cm“^ . This provides 
the potential to allow in situ creation of 1 km^ reflective 
surfaces with masses of 20 to 100 kg. For these, struc¬ 
ture and propulsion/control systems will likely be the 
dominant mass. Each reflector could be pointed based 
upon orbital position, and controlled independently and 
autonomously. 

We will not attempt a detailed discussion of the dy¬ 
namics of such an installation in this paper. However, 


we can utilize the planet’s L2 Lagrangian point to esti¬ 
mate the largest stable orbits in the effective planetary 
gravitational well. For a planet in orbit at the inner edge 
of t he HZ of an M8 star (as estimated below in Section 
|3.1| ), the Lagrange point is beyond 10 planetary radii, 
and is even more distant for earlier stellar types. There¬ 
fore, for mirror flotillas with orbits of modest extent (< 
5 planetary radii), the mirror orbits should be relatively 
stable and controllable. 

Providing a reflective area equivalent to Earth’s disk 
area (1.3x10® km^ ) requires deployment of more than 
130 million mirrors of 1 km^ each, or 1.3 million mirrors 
of 100 km^ . At 1000 kg km“^ , even a 100 km^ reflector 
lies within current human launch capabilities. However, 
even with optimistic ($2200/kg) estimates®, launch costs 
for such a fleet exceed $ 2 . 8 x 10 ®® USD 2015 . 

There will undoubtedly be significant environmental 
impacts due to pollution from manufacturing, launch ex¬ 
haust products and altering the energy balance of the 
planet. However, we leave the details of engineering a 
climatic balance, combating pollution, and the possible 
sacrifice of a millennium of economic output to the in¬ 
habitants of these worlds. 

It is worth highlighting certain considerations affecting 
implementation of dark-side illumination. For example, 
heating the dark side would lower the day-night tem¬ 
perature differential and slow day-to-night heat transfer, 
warming the day side. The mirrors would also obscure 
some starlight for day-side inhabitants. Optimally the 
mirrors could be turned edge-on to our transit line-of- 
sight quickly as they exit the zone of usefulness, in order 
to minimize obscuration, unnecessary exposure to solar 
radiation, and forces that would alter the orbit. How¬ 
ever, some obscuration may be desirable to reduce the 
stellar flux to the day side and keep planetary tempera¬ 
tures from rising excessively. Obscuration of the day side 
isn’t overly large even if the mirrors are not directed edge 
on, lying between ~ i?®/i?^ and ~ R^/Rf^ depending on 
the orbital arrangement, where Rm is the orbital radius 
of the mirror fleet. Satellites in a 3Rp orbit would de¬ 
crease the flux by approximately 4 to 11%, while a 10i?p 
orbit reduces flux by 0.1% to 1%. 

As a result of these and other considerations, the de¬ 
tailed distribution and orientation of mirrors within the 
orbiting fleet are likely quite complicated. The orbital 
distribution and variable incidence angle of the mirrors 
would affect the actual mirror surface area required, as 
well as the observed transmittance of the mirror fleet. 
However, for most of this paper we restrict ourselves to 
the simplest case and consider the observable impact of 
a constant-transmittance fleet of mirrors surrounding an 
extrasolar planet. Near the end of Section]^ we briefly 
consider alternative transmittance models, while defer¬ 
ring more detailed modeling of mirror fleet transmittance 
to another time. 

3. MODELING EFFECTS OF A CLOUD OF MIRRORS 

An orbiting mirror fleet partially blocks starlight just 
before, during, and after planetary transit. We initially 

® In 2015, launch costs to geosynchronous transfer or¬ 
bit on Falcon Heavy in 2015 were estimated at $14100/kg: 
http://www.spacex.com/about/capabilities, retrieved on May 16, 
2015 
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assume the fleet appears in cross-section as a constant- 
transparency ring around the planet. Although unlikely 
to be completely accurate, these models highlight the 
first-order effects of orbiting mirrors on transit light 
curves. We also neglect the increased light reflected from 
the now-illuminated planetary dark side, an assumption 
validated in Section [Ql 

We modeled fleets of orbiting mirrors in two 
ways. For our standard r uns, we adapted the EXO- 
FAST_OCCULTQUAD r outines dEastman et al .|20f3 ) based 
on the algorithms of iVlandel &: Agol 1 2(J(J2|)' i'his pack¬ 
age calculates light curves as a function of the normalized 
separation of the star and planet centers (z = d/i?star), 
which for circular orbits is related to time (t; 0 at transit 
center) by 


z = 


ap 

Rstai 


\/ (sinwt )2 -p (cosicoswt)^ 


where i is the orbit inclination and ap and w describe 
the radius and angular frequency of the planetary orbit. 

The EXOFAST routine is used to produce a light curve 
for a transit by a planet without mirrors (LGp), and 
again to produce a light curve for a transiting object the 
size of the mirror fleet {LCq)- The fraction of starlight 
blocked in each case is Bp = 1 — LCp and Bq = 1 — 
LCo, respectively. The light blocked by the transparent 
mirror ring (absorptance = a = 1 - transmittance) is 
therefore Bm = ct{Bo — Bp) and the light curve for a 
planet encircled by the mirror fleet is 

LCtot = l-BM-Bp = LCp{l-a) + aLCo ( 1 ) 

The absorptance a and optical depth r are nearly equiv¬ 
alent for a < 0 . 1 ; however this is not necessarily valid 
for our simulations. 

We assume the effective cross-sectional area of the 
constant-transparency mirror ring (aAm) equals the 
cross-sectional area of the planet (Ap), resulting in an 
absorptance 


a = 


Ap 

Am, 


ttB-p 


1 


Tr{R‘^ - R%) 




( 2 ) 


This is reasonable because any mirrors between Earth 
and the planet are not illuminated, so do not contribute 
to the area needed to illuminate the dark side at day- 
side levels. We also assume mirrors are turned edge-on 
quickly as they exit the zone of usefulness, i.e. Any 
mirrors in a position that renders them incapable of il¬ 
luminating the planet’s dark side are edge-on so don’t 
contribute to the fleet’s absorptance. 

Although our standard simulations assume the mir¬ 
ror fleet provides constant absorptance, we wished to 
explore the importance of this simplification. We there¬ 
fore created software to generate a two-dimensional array 
representing the absorptance of the planet and mirror 
fleet, and another representing the star. For all runs, 
the planet radius was set to 100 pixels. We then shifted 
the relative positions of the two arrays, calculated the 
transmitted light for each location, and interpolated the 


transit depth, within the range expected due to pixeliza- 
tion errors. We used this approach to simulate transits 
for configur ation s with non-constant absorptance profiles 
(see Section 4.3). 

3.1. Parameterization 

Table [l] summarizes the properties of stars used in our 
transit models. Columns 1-5 co ntain the spectra l type 


and basic stellar properties from Zombeck (1990). Val¬ 
ues for surface gravity log(g) were estimated two ways. 
For Method A we retrieved information on Kepler stars 
hosting exoplanets from the Exoplanet Data Explorer 


(Wright et al. 2011) at http://exoplanets.org/table, fit 
the data with log{g) = Ci+C 2 leS, then utilized this func¬ 
tion to estimate log{g). For Method B we scaled the sur¬ 
face gravity from the solar value, assuming g cx M/R?. 
In several test runs, the differences in the resulting light 
curves were at most 1 % of the transit depth, and gener¬ 
ally less than 0.2%. For all of our standard runs we used 
the values from Method A, as shown in column 6 . 

For our standard runs, we assumed heavy-element 
abundances half that of the Sun (log([Fe/H]ratio) = 
—0.3), and calculated quadratic limb-darkening pa- 
rameters 7 i and 72 using the interface of |East- | 
|man et al.| (|2013|) ( http://astroutils.astronomy.ohio- 
state.edu/exofast/l imbdark.shtml), which i nterpolates 


over the values of Claret & Bloemen (2011). The rel¬ 
atively low element abundance value was chosen because 
we anticipate advanced civilizations around older stars 
which tend to be metal poor. 

We assumed circular planetary orbits in all our mod¬ 
els. Edge-on orbits were standard, but we also explored 
other orbit inclinations. To determine appropriate orbit 
sizes, we estimated each star’s habitab l e zone using the 
optimistic values of [Kopparapu et al.] ( |2013| ) , accessed 
via http: / / depts.washmgton.edu / naivpl / content/hz- 
calculator. We note that substantial disagreement exists 
on the precise range of orbits which lie in a star’s HZ . 


who compare the conservative (rather than optimistic) 

choices made by Kopparapu et al. ( 

20131) with those of 

other authors (Kastmg et al. 


Zsom et al. 20131 

Pierrehumbert Gaidos 2011 

). However, in this paper 


we are concerned mainly with illustrating the potential 
impacts of a fleet of orbiting mirrors on transit light 
curves, using one plausible choice for the HZ. 

Table summarizes the properties of planets used in 
our transit models. For each stellar spectral type (col¬ 
umn 1 ), the inner edge (HZin) and center (HZ^id) of the 
chosen HZ are indicated in columns 2 and 5. Correspond¬ 
ing planetary orbit periods for circular orbits (Porb.in 
and Porb.mid), determined from stellar mass and orbit 
size, are in columns 3 and 6 . Because we anticipate mir¬ 
ror fleets might be implemented around planets in syn¬ 
chronous rotation, columns 4 and 7 contain estimates of 
the tidal synchronization timescale (Tgyn) for Earth-like 
(M = Me, R = Re) planets in these orbits. Plan¬ 
ets which would not synchronize within ~ 10 Gyr are 
marked with an asterisk (*). 


with EXOFAST. This brute-force modeling approach is 

Leconte 

(|2013|) 

‘syn wc L 

based on 

significantly slower than that already described, but in 

1 Darwin 

1880 f 

Hut 

1981 


based on the equilibrium theo ry ot tides 


ated by the two methods differed by at most 0.05% of the cording to Lunar Laser Ranging experiments of Neron de 











































































Technology in Light Curves 


5 


Table 1 

Stellar Parameters 


Sp Type 

TMK) 

M/M,u„ 

R/Rsun 

L/Lgun 

log(3) 

M8 

2660 

0.10 

0.13 

0.0008 

4.96 

Ms 

3120 

0.21 

0.32 

0.0079 

4.87 

Mo 

3920 

0.47 

0.63 

0.063 

4.72 

Ks 

4410 

0.69 

0.74 

0.16 

4.62 

KO 

5240 

0.78 

0.85 

0.40 

4.46 

CS 

5610 

0.93 

0.93 

0.79 

4.39 

C2 

5780 

1.00 

1.00 

1.00 

4.35 

Co 

5920 

1.10 

1.05 

1.26 

4.33 


Table 2 

Planet Parameters 


Sp Type 

HZi„ (AU) 

Porb.in (d) 

-Tayn.in (Gyr) 

HZ,„id (AU) 

Porb.mid (d) 

‘^svn.mid (Gyr) 

M8 

0.023 

3.27 

3 xlO-'^ 

0.043 

9.24 

1 X 10"® 

Ms 

0.073 

14.8 

7 xlO-^ 

0.13 

37.4 

2 xlO-'^ 

Mo 

0.20 

47.7 

6 xlO”'^ 

0.36 

115 

0.2 

K5 

0.32 

79.6 

5 xlO-2 

0.55 

179 

1 

KO 

0.48 

138 

0.4 

0.83 

313 

ll' 

G5 

0.67 

208 

2 

1.13 

455 

50 ' 

G2 

0.75 

237 

4 

1.26 

517 

00' 

Go 

0.84 

268 

6 

1.40 

577 

100' 


Note. — * Earth-like planets in these orbits do not synchronize their rotation within 10 Gyr. 


Surgy & Laskar (1997), which yields 


r,y„(Gyr) = 21.71 


(Mp/Mi5)a(AU)6 


(3) 


Terrestrial planets at the inner edge of the chosen habit¬ 
able zone (HZin) synchronize their rotation within 6 Gyr 
for all G, K and M main sequence stars. However, Tsyn 
is longer for terrestrial planets elsewhere in the HZ, For 
example, for planets at its center (HZmid), synchronous 
rotation is achieved within ~ 1 Gyr for main sequence 
stars K5 and later, but not within 10 Gyr for stars KO and 
earlier. Tidal synchronizat ion time-scales are u ncertain 
at best. Estimates made by Kasting et al. (1993) suggest 
that for planets orbiting at HZ^id, only tnose around M 
stars would be in synchronous rotation at the current 
epoch, while those orbiting at HZin around most K stars 
could also be in this state. Although detailed estimates 
vary, planets in an M star’s HZ are extremely likely to 
be in synchronous rotation after a several Gyr biological 
evolution timescale. However note that a planet’s climate 
and habitability may also be altered by effects such as 
tidal heating or orbital migration (see Barnes et al.|200^ 


2013). As they are most likely to be in synchronous ro- 


tation, we focus on planets around M stars for the ma¬ 
jority of our models, but do explore other options, since 
Earth-like planets orbiting a G star near HZjn might be in 
synchronous rotation after sufficient time, and as noted 
earlier, civilizations on Earth analogues might choose to 
utilize a fleet of mirrors for reasons other than dark-side 
illumination of a synchronously rotating planet. 

In each orbit, we simulated transits for planetary radii 
of i?p/i?Earth = [0.5, 1, 2], Corresponding to M/Me = 
[0.125,1,8] if all have the same average density as Earth, 
or M/Me = [0. 53,1, 5.1] using the relationships offWeiss 
& Marcy (2014). Effects on the transit due to brightness 
ot the planet’s dark side were ignored. For each planet. 


we simulated mirror fleets extending to three different 
distances from the planet center: Rm/Rp = [2, 3, 10]. 
Mirror fleets that are substantially smaller {RMRp = 
•\/2) will be nearly opaque according to equation^ while 
those that are substantially larger may have less stable 
orbits and may suffer from inefficiencies in redirecting 
starlight. 

For all runs, in addition to determining the light curve 
for the planet surrounded by its mirror fleet (P + M), we 
also determined the light curve for the planet alone (P) 
and for a larger planet (PP) sans mirrors that produces 
the same eclipse depth as P + M. 

4. RESULTS OF SIMULATIONS 

Overall, mirror-induced perturbations to the transit 
shape are similar for all situations. We begin by pre¬ 
senting one model as a reference, then explore the effects 
of varying stellar, planetary, and mirror fleet parameters. 
We chose a large terrestrial planet orbiting an M5 star 
as our reference, because these stars and planetary sizes 
are relatively common, and a relatively large planet will 
create a deeper transit with a longer entrance time, mak¬ 
ing it easier to analyze the transition region potentially 
affected by the mirror fleet. 

The top image in Figure [^displays our reference sim¬ 
ulation: a planet with a radius Rp = 2PEarth, sur¬ 
rounded by a constant-absorptance fleet of mirrors with 
Rm = 3Pp, transiting in front of an M5 star. The solid 
line in the left panel of Figure [^displays the simulated 
light curve for the planet and its neet of mirrors (P+M), 
with the dashed line representing the planet-only transit 
curve (P). A transit of the same depth as P + M but 
caused by an isolated larger planet (PP, as in the bottom 
of Figure 0 is shown by the dotted line. 

Gontaimng many individual reflectors, the translucent 
mirror fleet blocks some starlight (see Figure and 
deepens the transit. This P + M case can be aistin- 
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Figure 2. Top panel. Simulated image of a planet with Rp = 
2/?Earth? surrounded by a constant-absorptance mirror fleet with 
R^ = SRp, transiting in front of an M5 star. Only a strip of the 
right half of the star is shown. Bottom panel Same as above, but 
for an isolated planet large enough to produce a transit of the same 
depth. 



z z 


Figure 3. Left panel Transit light curves for a planet with Rp = 
2/?Earth Orbiting an M5 star (dashed line: P), the same planet 
surrounded by a constant-absorptance mirror fleet extending to 3 
planetary radii (solid line: P -h M), and a planet without mirrors 
large enough to produce a light curve with the same depth of transit 
(dotted line: LP). Right panel Difference between the mirror fleet 
transit light curve {P-\-M) and the one for the solitary large planet 
(LP), relative to the stellar intensity. 

guished from that of an isolated larger planet (LP) in 
the timing and shape of the transit ingress and egress 
light curve. In the initial phases of the transit, only 
the mirrors block starlight, so the decrease occurs rel¬ 
atively gradually. When the planet itself enters transit, 
the light curve decreases more rapidly. Later in the se¬ 
quence, the dimming slows with the planet fully in front 
of the star while the semi-transparent halo of mirrors con¬ 
tinues to occlude more starlight. These effects are seen 
in the left-hand panel of Figurej^ although the signature 
is complicated by the presence of stellar limb-darkening. 
Note that once the mirrors are fully in front of the star, 
as depicted in Figure!^ the transit light curve becomes 
very similar to that ofa larger isolated planet, with the 
small differences caused by variations in stellar brightness 
across the differing areas of occlusion. The right panel 
shows the difference of the light curves for the planet plus 



Figure 4. Impact on light curve of abundance-induced limb- 
darkening variations compared with that of a mirror fleet. The 
solid line displays the difference (relative to stellar intensity) be¬ 
tween light curves produced by a planet surrounded by a mirror 
fleet and an isolated large planet, both transiting an identical M5 
star with half-solar abundance. The other two lines display the 
difference between the isolated large planet light curve for half¬ 
solar abundance and the light curves for isolated planets producing 
transits of the same depth for tenth-solar abundance (dotted) and 
full-solar abundance (dashed). 

mirrors transit and the larger planet transit, highlighting 
the effects as the planet enters and exits its transit. 

Although a fleet of mirrors orbiting an exoplanet af¬ 
fects both the time-span and depth of a transit, it would 
not affect the timing or depth of the planetary eclipse 
optical light curve. Once again, a difference in transit 
and eclipse time-span could suggest the presence of mir¬ 
rors to illuminate a planetary dark side, but the same 
caveat regarding elliptical planetary orbits applies here 
as for the single large mirror. We assume the mirror fleet 
is designed for maximum efficiency, so all starlight is di¬ 
rected onto the planet and none is re-directed towards 
us during the eclipse. However, if the fleet does not in¬ 
clude efficiency as a primary design constraint, we might 
observe a temporary increase in brightness just before or 
after the eclipse phase due to spillover effects. 

We anticipate no detectable effects due to thermal 
emission from either a single large mirror or an armada 
of orbiting mirrors. We assume the back sides are cold, 
as mirror reflection efficiency is high and no bright IR 
sources illuminate them from behind. Since light mir¬ 
rors are very thin, the front and back are at nearly the 
same cold temperature. 

4.1. Effects of Varying Stellar Abundance or Spectral 

Type 

Since the mirror fleet mainly affects transit entry and 
exit, we questioned whether differences in limb-darkening 
could exhibit the same signature. One source of un¬ 
certainty in limb-darkening arises from abundance vari¬ 
ations. In runs for identical planets (Rp = 2i?Earth, 
Rm = iRp), varying the stellar abundance from half¬ 
solar to tenth-solar or to full-solar impacted the light 
curve depth by 0.5% to 4%, depending on stellar type. 
This is of similar magnitude as the effects of a mirror 
fleet. We therefore created light curves for isolated plan¬ 
ets transiting M5 stars with tenth-solar and full-solar 
abundances, with sizes generating transit depths equal to 
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that for the reference simulation. The solid line in Figure 
15 shows the difference of the mirror plus planet and large 
^anet light curves from Figure [5 while the dotted and 
dashed lines show the difference m the light curves for the 
isolated planets for various stellar abundances. Different 
stellar abundances cannot be mistaken for the effects of 
a transiting planet surrounded by a fleet of mirrors, as 
limb-darkening variations have a fundamentally different 
signature. 

Although not displayed, we simulated transits of plan¬ 
ets orbiting the full range of stellar types displayed in 
Tabled Transit light curve details differ somewhat for 
stars 01 earlier spectral type. In particular: (1) Differ¬ 
ences in limb-darkening affect the detailed shape of the 
light curve near entry and exit. (2) For larger stars, 
the same planet plus mirror fleet blocks a smaller frac¬ 
tion of light, resulting in a shallower transit, and the 
mirror fleet’s effects last a smaller fraction of the tran¬ 
sit time-span. Compared with an isolated planet transit 
of the same depth, the impact of the mirrors is of order 
10“"^/star for an MO star, roughly a factor of 2 smaller for 
a G2 star, and several times larger for an M5 star. (3) 
Longer orbital periods, corresponding to the larger HZ 
orbits around earlier-type stars, increase the time span 
of transits for these planets. Despite these differences, 
simulating transits of other stars does not significantly 
change the bipolar signature of the mirror fleet in transit 
residuals. In addition, the maximum residuals remain 
approximately the same relative to transit depth for all 
stellar types (e.g. ^4-6% for Rp = 2/?Earthj Rm = SRp)- 

4.2. Varying Mirror Fleet Extent, Planetary Size or 
Orbit Inclination 

The top panel of Figure displays light curves for 
transits that are identical to the reference simulation 
apart from the mirror fleet’s extent. Larger fleets affect 
the light curve for a larger fraction of the transit dura¬ 
tion, with more dramatic departures from the isolated 
planet situation. The dashed curve [Rm = 10i?p) has a 
slightly different transit depth because Rp/Rstar = 0.057 
so Rm/Rstar = 0.57; the mirror fleet extends across a sig¬ 
nificant fraction of the stellar surface, which has differ¬ 
ent intensities at different locations due limb-darkening. 
Because all mirror fleets have the same effective area, 
the larger fleet has a lower absorptance. So at z = 0 
the smaller, more opaque systems block light primarily 
from the star’s brighter central regions, while the largest, 
most transparent one transmits more of that light, in¬ 
stead blocking less intense outer portions of the star. The 
bottom panel illustrates that, all else being equal, mirror 
fleets with large spatial extent would be more apparent in 
transit light curve abnormalities than those with smaller 
spatial extent, both temporally and in the magnitude of 
the effect on transit residuals. 

For mirror fleets as in the reference model, we explored 
the effects of varying planet size. Figure [5 displays the 
the difference between the P + M and LP light curves, 
relative to transit depth /star — /min) for three sizes of 
terrestrial planets, each orbiting at HZiajd around an M5 
star. Smaller planets result in transits of substantially 
shallower depth, however in all cases the mirror fleet has 
a significant impact (a few percent of the transit depth) 
on the entrance and exit light curves. The duration of 




Figure 5. Top panel: Transit light curves that result when a 
planet with Rp = 2/?Earth? located in the middle of the star’s HZ, 
passes in front of an M5 star. In all cases the planet is surrounded 
by a constant-absorptance mirror fleet, with Rm = SRp (solid), 
Rj^ = 2Rp (dotted), or Rm = 10i?p (dashed). Bottom panel: 
Difference between the mirror fleet transit light curve (P + M) and 
the one for a solitary larger planet (LP) that would produce the 
same depth of transit, relative to the stellar intensity, for the same 
situations. 



t (days) 


Figure 6. Impact of a fleet of mirrors with Rm = SRp on transits 
for planets of various sizes orbiting in the middle of an M5 star’s 
HZ. The lines show the difference between the light curve result¬ 
ing from the planet plus mirror (P -h M) transit and one of the 
same depth resulting from an isolated larger planet {LP) transit 
for planets with Pp/PEarth = 0-5 (dashed), 1.0 (dotted), and 2.0 
(solid). 
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Figure 7. Top panel: Transit light curves that result when a 
planet with Rp = 2/?Earth? located in the middle of the star’s HZ, 
passes in front of an M5 star. In all cases the planet is surrounded 
by a constant-absorptance mirror fleet with Rm = lORp. The 
four curves result from different values for the impact parameter: 
•^min = 0 (solid), 2:min = 0.75 (dashed), = 1 (dotted), and 

= 1-25 (dash-dot). Bottom panel: Difference between the 
mirror fleet transit light curve (P -h M) and the one for a solitary 
larger planet (LP) that would produce the same depth of transit, 
relative to the stellar intensity, for the same situations. 

largest residuals is shorter for smaller planets. 

All previous results are for edge-on orbits, so the mini¬ 
mum projected distance between star and planet centers 
is Zmin = 0. The top panel of Figure [^compares this sit¬ 
uation (solid line) with orbits that are not exactly edge- 
on (zmin 7 ^ 0). For these, a planet with Rp = 2 Re and 
Rm = 10i?p is located at HZ^id around an M5 star. The 
dashed line planet transits at Zmin = 0.75, so part of the 
mirror fleet never occludes the star {Rm/Rstai = 0.57). 
In the dotted line case, the planet center crosses the edge 
of the star (zmin = 1), and for the dash-dot light curve 
(zmin = 1.25) the planet itself never crosses in front of the 
star, but the mirrors still create a transit. For each im¬ 
pact parameter Zmin, the bottom panel of Figureshows 
the difference of the P + M and LP light curves. This 
figure reveals that as Zmin increases, the mirror effects 
are evident during a greater fraction of the transit time- 
scale. As long as more than half of the planet crosses 
in front of the star, i.e. z^in < 1; the transit residuals 
display a bipolar signature of the mirror fleet’s presence. 
For Zmin > 1, the residuals are no longer bipolar, but 
the mirror fleet produces transits that last much longer 
than would be appropriate for the transit depth in the 


absence of mirrors, and the light curve minimum is not 
at transit center. 


4.3. Effects of Different Mirror Density Distributions 

Having explored the first-order effects of a fleet of 
mirrors, we considered alternative transmittance models 
based on actual distributions of steerable mirrors, which 
required our brute-force approach. The simplest such sit¬ 
uation is with mirrors evenly distributed (density Pmirr) 
within a spherical shell (i?^ < r < R^) around the planet. 
Although not necessarily realistic, this scenario lets us 
explore the importance of our constant-absorptance as¬ 
sumption. We calculated the absorptance at each pro¬ 
jected radial distance based on the projected density of 
the fleet: 

Pproj = 2pmirr( V R^ - r^ - Rf - r^) {r <Ri) 

Pproj = 2/9mirr(V^o “ {Ri < r < Ro) 

The absorptance was again normalized so the total light 
intercepted equals that normally incident on the planet, 
while simultaneously ensuring absorptance a < 1. We 
considered two extremes: a thick uniform shell of mir¬ 
rors extending outward from just above the planet’s sur¬ 
face, and mirrors all at nearly the same altitude. In the 
former case, the mirror fleet’s collective absorptance is 
largest near the planet’s surface (near Ri), gradually de¬ 
creasing to zero at the fleet’s outer extent. For the thin- 
shell case, the absorptance is small near the planet’s sur¬ 
face, increasing with distance, then sharply peaked near 
the fleet’s altitude. Figure [^illustrates the extreme dif¬ 
ferences in absorptance profiles for the thick-shell (top 
panel) and thin-shell (bottom panel) cases. 

These simulations illustrate the relative insensitivity 
of transit light curves to varying mirror-fleet absorp¬ 
tance profiles. Figure [^ compares transit light curves 
for these two cases with previously discussed results. 
The constant-absorptance and large planet light curves 
are the same as those shown in Figure The dotted 
fine results from a thin shell of mirrors (i?^ = 2.99i?p, 
Rm = 3i?p), and the dashed line from a thick shell of mir¬ 
rors extending from Ri = 1.01i?p to Rm = 3i?p. The in¬ 
ner radii in both extremes were chosen based on the reso¬ 
lution of our brute-force simulations (planet radius = 100 
pixels). From this figure, it is apparent that light curves 
for various absorptance distributions are more similar to 
each other than they are to that for an isolated planet. 
Our simple constant-absorptance model therefore cap¬ 
tures the dominant effects of this method of dark-side 
illumination. 

5. PROSPECTS FOR DETECTION 

We now consider the potential detectability of such 
mirror installations, focusing on the predicted {P -I- M) 
light curve for an Earth-sized planet {Rp = i?Earth) or¬ 
biting a G2 star at HZin, surrounded by a large {Rm = 
lORp) constant-absorptance fleet of mirrors. We chose 
this situation in part because the only known civiliza¬ 
tion evolved on such a planet. Older M stars are very 
common, and their potentially habitable terrestrial plan¬ 
ets are likely in synchronous rotation. However they are 
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Figure 8. Top panel. Simulated image for the same situation as 
Figure but for a planet surrounded by a constant-density fleet 
of mirrors in a thick shell extending from just above the planet’s 
surface (from 1.01i?p to 3Rp). Bottom panel. Same as above, but 
the mirrors are distributed in a thin shell extending from 2.99Rp 
to 3Rp. 



Figure 9. Transit light curves for different absorptances. The 
solid line displays the light curve of a transit produced by a planet 
with Rp = 2i?Earthi s-nd a constant-absorptance fleet of mirrors 
extending outward to Rm = 3Rp, crossing in front of the center of 
an M5 star with half-solar abundance. The dotted and dashed lines 
display the transit light curve produced by thick {Ri = l.OlfJp) and 
thin {Ri = 2.99Rp) constant-density shells of mirrors, respectively, 
each with an outer radius of 3Rp. The dash-dot line is the light 
curve for an isolated planet large enough to produce a transit of 
the same depth. 


fainter, so not necessarily better candidates when seeking 
signatures of planetary-scale engineering. Recall that a 
planet near HZin for a G2 star m ight be synchronously 
rotating at this time (see Section 3.11, or that civiliza¬ 
tions might choose to affect the climate of bodies other 
than their planet of origin. 


5.1. Detectability with Kepler 

Until now we have neglected the reality that a mirror 
fleet’s effects will be partially suppressed by transit fit¬ 
ting routines through variations in orbit inclination or 
stellar limb-darkening. To assess this, we re-sampled our 
P + M model onto Kepler's short cadence (SC) timescale 
for 3 transits. We used the Transit Analysis Package 
(TAP, Gazak et al. (2012)) to fit an isolated planet model 


to the simulated light curve, restricting limb-darkening 
parameters to a range appropriate for stars with effec¬ 
tive temperatures within 750K of a G2 star. Because the 
fleet is translucent (a ^ 0.01 for Rm = lORp), predicted 
residuals (~ 10“®/star) are substantially small er than 
those for the Jupiter-scale opaque structures of [Arnold] 
(20051, and thus too small to detect in Kepler data. 

id confirm this, we simulated “best-case scenario” Ke¬ 
pler data for our P -\- M model and compared the resid¬ 
uals for the previously determined isolated planet model 
with the difference between the P M model and the 
isolated-planet TAP model determined above (hereafter 
TAP IPM). Given the relatively long orbital period (^ 
237d) for a planet orbiting a G2 star at HZi„ (a = 0.75 
AU), at most 6 transits of data are available. Kepler 10 
is a relatively bright (Kp = 10.96) G2 star in the Kepler 
sample, with instrumental and photon noise estimated 
by the Kepler pipeline as ~200 parts-per-million (ppm) 
for short-cadence data. We simulated 6 transits of SC 
data with Gaussian white noise at this n oise level, al¬ 
though we note that Batalha et al. (2011) report more 
long-cadence noise (^62 ppm) than the Kepler pipeline 
estimate (^ 36 ppm). 

The top panel of Figure [T0| displays the phase-folded 
simulated data, where we retain the folded time scale for 
convenience in later analysis and t=0 represents the cen¬ 
ter of each transit. For clarity, only half of the transit is 
displayed. The superimposed solid line is the TAP IPM 
described above, i.e. the one that best fits the noiseless 
P -\- M model used to simulated the data. The residu¬ 
als between our noisy data and the TAP IPM, smoothed 
for visual clarity, are shown in the middle panel. They 
are overlaid (solid line) by the smoothed difference be¬ 
tween our P-\-M model and the TAP IPM. Clearly, noise 
in these data overwhelms any potential mirror fleet sig¬ 
nal. In the bottom panel, the convolution of the (un¬ 
smoothed) simulated and predicted residuals confirms 
they are uncorrelated. 

M stars are substantially dimmer at visible wave¬ 
lengths, resulting in noisier data. A typical Kepler cata¬ 
log star has a magnitude of ~ 14, based on i nformation 


retrieved from the Exoplanet Data Explorer (Wright et 
al.|2011| , but the median M star is fainter than Kp = 15. 


Kepler 138 , hosting multiple terrestrial planets (Kipping 
2014 Rowe et al.|2014 |, is an unusually bright A e- 


et al. 


pier iVIU-dwart at Kp = 12.93, but still two magnitudes 
fainter than the G2 star simulated below. Thus the ex¬ 
pected noise level would be ^ 2.5 times greater. Be¬ 
cause of the shorter orbital period for a planet at HZin 
around an MO star, up to 5 times more transits of Kepler 
data could be available, mitigating the additional noise. 
However mirror fleets around these stars would still be 
undetectable in Kepler data. 


5.2. Future Prospects for Detection 
The James Webb Space Telescope {JWST) will have 
a collecting area of 25 m^, 35 times larger than Kepler 
(0.708 m^). Using area, quantum efficiency, and filter 
tran s mission parameters provided in|Van Cleve & Cald-1 
|well| (|2009D and http;// ircamera.as.arizona.edu/ 
nircam/features.html, along with the Planck black- 
body function for a G2 star, we computed the expected 
count rate of the NIRCAM F150W2 filter relative to the 
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Figure 10. Top panel: Simulated phase-folded data with 200 
ppm noise and Kepler SC sampling for 6 transits of a planet with 
Rp = -REarth Rm = lORp, Orbiting around a G2 star at 0.75 
AU. Only half of the transit is displayed. The solid line is the best- 
fit isolated planet model produced by TAP analysis of the noiseless 
P + M model. Middle panel: Smoothed simulated residuals (dots) 
of the simulated data and TAP model shown in the top panel, and 
predicted residuals (solid line) based on the difference between the 
isolated planet TAP model and the simulated P M model. The 
smoothing length for both was determined by half of the time- 
scale (defined in the text). Bottom panel: Convolution of the 
unsmoothed simulated and predicted residuals. The vertical dotted 
lines indicate ±4^/2 > relate to the statistical analysis described 
in the text regarding Figure m 

Kepler full bandpass. The JWST F150W2 count rate 
will be ^ 45 times that for Kepler, reducing noise levels 
to ~ 200/6.7 = 30 ppm. 

M stars emit a larger fraction of their light at infrared 
wavelengths, so we used the effective temperatures and 
radii of Table to calculate the relative expected JWST 
photon rates for G2, MO, and M5 stars at equal distances. 
In the F150W2 bandpass, we expect 6.7 and 53 times 
more photons from a G2 star than for MO and M5 stars, 
respectively. Data for M stars at similar distances will 
be noisier, requiring more transits to achieve detections, 
so only the closest and brightest M stars will be feasible 
targets for JWST. 

As before, we simulated data with short-cadence sam¬ 
pling for an Earth-like planet orbiting at a G2 star’s HZin 
using our P + M model, but at JWSTs expected noise 
level of 30 ppm. We compared the residuals between the 
simulated data and the TAP IPM with the difference be¬ 
tween the P + M model (used to simulate the data) and 
the TAP IPM. We did this comparison for 100 separate 




Figure 11. Same as Figure p^ but for 4 transits of simulated data 
with a noise level of 30 ppm, as anticipated for JWST. 

realizations of noise for each of 2,3,4,5, and 10 transits of 
JWST data. Although we would ideally base our com¬ 
parisons on isolated-planet TAP models fit to each simu¬ 
lated noise realization, this simplified first-order analysis 
is sufficient to illustrate the key points regarding mirror 
fleet detectability. 

The top panel of Figure [TT] shows the phase-folded data 
for a 4-transit noise realization, along with the TAP IPM. 
For visual clarity, only half of the transit is displayed; 
however all analyses were performed using the full-transit 
simulated data. The smoothed residuals in the middle 
panel reveal the signature of the mirror fleet; residuals of 
the simulated data and TAP IPM (dots) are compared 
with predicted residuals (solid line; difference between 
simulated noiseless P + M model and TAP IPM). Both 
are smoothed by half of the time-scale ti /2 (defined be¬ 
low) for visual clarity. The convolution in the bottom 
panel demonstrates the significant correlation between 
the (unsmoothed) simulated residuals and the predicted 
model difference. Vertical dotted lines delineate points 
satisfying |At| < ti/ 2 , where < 1/2 is the time scale for 
the isolated planet model light curve to decrease from 
maximum to half of the transit depth. The peak is 5.49(7 
above the convolution function’s mean, with both mean 
and RMS {a) determined by data outside the region de¬ 
lineated by 1 ^ 1 / 2 !, i.e. where we do not expect a signifi¬ 
cant signal. We display this instance of noise because its 
convolution peak significance (5.49 a) equals the median 
significance for the 100 separate 4-transit noise realiza¬ 
tions. 
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In the 165 points with |At| < t-iji, the probability of 
the convolution having a peak of at least this significance 
occurring by chance is 4.4 x 10“^%. Of the 100 instances 
of 4 transits with white Gaussian noise of 30 ppm, the 
convolution peak was at least 4cr 94 times. In our 5- 
transit simulations, 99/100 have convolution peaks 4 (t 
or higher, with a median of 6.10cr, and a probability of 
random occurrence no greater than ^ 2% in all 100 re¬ 
alizations. Thus for a star like Kepler-10, we anticipate 
successful detection of such a mirror fleet with data for 
only a moderate number of transits. 

Could the effects of dark-side illumination be de¬ 
tectable with fewer transits? For the 100 separate 3- 
transit realizations, the median significance of the con¬ 
volution peak was 4.78cr, with 83/100 occurring at Aa or 
above. For the 2-transit realizations, the median peak 
significance was 3.90tT, with only 44/100 reaching the 4cr 
threshold. We conclude that at this 30-ppm noise level, 
spurious low-significance detections become increasingly 
likely with fewer than 4 transits of data. 

An initial impression of the 4-transit residuals as shown 
might suggest the isolated planet fit is adequate, so mir¬ 
ror fleet signatures could go unnoticed without a correla¬ 
tion analysis similar to that carried out above. However, 
for 10-transit simulations, the convolution analysis re¬ 
sults in a peak above 4cr 100% of the time, with a proba¬ 
bility of the weakest peak occurring by chance of <0.06%. 
Figure [T^ displays the results of a 10-transit simulation 
with a peak of median significance (7.74cr). In this situ¬ 
ation, the simulated smoothed residuals do indicate that 
the isolated planet model is insufficient, and the high sig¬ 
nificance of this signature is made apparent by the strong 
peak in the convolution shown in the bottom panel. The 
probability of this peak occurring at random is only 2.7 
xl0-i°%. 

Residuals could be further suppressed if the isolated- 
planet TAP models do not restrict the range of limb- 
darkening parameters to values appropriate for similar 
stars. However, in the majority of cases the stellar tem¬ 
perature and spectral type are known sufficiently well to 
support the use of some restrictions. 

The Transiting Exoplanet Survey Satellite ( TESS), ob¬ 
serving at 600-1000 nm, may identify transiting exoplan¬ 
ets for which targeted JWST observations would be suf¬ 
ficient for our purpose (Ricker 2014, Ricker et al. 2014). 
TESS will look at stars 30 to 100 times brighter than 
Kepler (5 < m < 12), across the full sky. Primar¬ 

ily sensitive to orbital periods up to 60 days, it can de¬ 
tect longer orbital periods in regions of the sky coin¬ 
ciding with JWSTs continuous viewing zone (CVZ). It 
will study primarily F, G, and K stars, but also all of 
the M dwarfs within 200 light years. It is anticipated 
that TESS will identify ^ 300 transiting super-Earths 
and tens of Earth-sized planets, with ~ 1/3 of them in 
or near the GVZ. TESS might well discover Earth-sized 
exoplanets in the habitable zones of K or M stars bright 
enough for detailed follow-up observations with JWST or 
future ground-based 30-m class telescopes. Since high- 
quality transit data on such planets are of great interest 
for studying properties such as their atmospheres, data 
adequate for the task of identifying or ruling out fleets 
of orbiting mirrors will become available. 




Figure 12. Same as Figure |10| but for 10 transits of simulated 
data with a noise level of 30 ppm, as anticipated for JWST. 


5.3. Second-Order Effects & Potential False Positives 

Since orbiting mirrors illuminate the planet’s dark side, 
we explored the importance of planetary reflection of this 
light. We assume the planetary contribution to photons 
received is constant, and that the amount of starlight 
reaching the day side equals that reaching the “dark” side 
via mirrors. The planet to stellar flux ratio is therefore 
of order 



The first factor represents the fraction of total starlight 
intercepted by the cross-sectional area of the planet. A 
is the planet’s albedo, i.e. the fraction of starlight that 
is scattered, with division by 2Tr because it scatters into 
27r sr. The final factor accounts for dimming of scat¬ 
tered planetary light by the mirror fleet; absorptance is 
reduced because only half the mirrors lie on our side of 
the planet. Greater effects occur for planets that are 
large relative to the stellar size. For small M8 stars, a 
planet with Rp = 2i?Earth has radius Rp — 0.14i?M8- 
For a mirror fleet extending to Rm = 3i?p, the absorp¬ 
tance (given our assumptions) is a = 0.125. Assuming 
a planetary albedo comparable to Earth’s (A = 0.3), 
this contribution alters the observed transit depth by at 
most 0.089%, significantly les s tha n the 4-6% due to mir¬ 
ror fleets (see end of Section 4.1 1 . Thus scattered light 
from the now-illuminated dark side produces only a small 
effect on a transit light curve. 
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Since JWST and other future telescopes will observe 
primarily in the infrared, we must consider thermal ra¬ 
diation from the now-warmed dark side. The planet to 
stellar flux ratio at any particular wavelength is 

Fp{\) ^ eB^jTp) / Rp y 

^star('^) '^A(^star) \star / 


where B\{T) is the Planck thermal blackbody radia¬ 
tion, and e is the planet’s emissivity. An upper limit 
is achieved by assuming e = 1 and choosing the coolest, 
smallest star (M8). For an M8 star, i?p/i?star = 0.14 
and Tstar = 2660K; we assume a habitable planet with 
Tp = 300K. In the JIFS'Tbandpass (0.6p,m-28.5/rm), the 
flux ratio reaches a maximum of 0.09% at a wavelength of 
28.5/im. At the wavelength of strongest planetary emis¬ 
sion (~ lO^m), the flux ratio is only 0.01%. The change 
in transit depth is significantly smaller than the 4-6% ef¬ 
fects of orbiting mirrors. Even if Tp = lOOOK, planetary 
thermal radiation only alters the transit depth by 0.7%. 

In general, any potential detection of intelligence must 
be confirmed using an entirely different method. We 
must consider whether natural phenomena could mimic 
signs of intelligence. In particular, a face-on ring system 
is extremely similar to the translucent projected annu¬ 
lus of our putative fleet of mirrors. Several gr oups have 
modeled transit detectability of ring systems dBarnes & 


Fortney 2004 Ohta et al. 2009 Tusnski & Vallb 2011), 


generally focusing on gas giant rings. JNo known terres¬ 
trial planet has a ring system, and we question the stabil¬ 
ity of such a system around a terrestrial planet i n tidally- 
ind uced s ynchronous rotation, particularly since |Ohta et 
ST] p009| found that rings around all close-in planets 
are likely to be short-lived and nearly edge-on, render¬ 
ing them more difficult to detect. More importantly, ring 
particles are likely warmer than the mirrors, emitting sig¬ 
nificantly in the IR thus reducing the IR transit depth. 
Comparison of IR and optical light curves will therefore 
help discriminate between natural and engineered sys¬ 
tems, given sufficient sensitivity. Ring systems may also 
be distinguished from mirror fleets by the effects of for¬ 
ward scattering of light by ring particles, causing flux to 


brighten just be fore and after planetary transit (Barnes 
fc Fortn^|2004 ). 


Planets m an M star’s HZ might experience at¬ 
mospheric erosion, due to strong stellar activity and 
weak planetary magnetic moments (e.g. Lammer et 
al. 2007 and discussions in Scalo et al. 2007 and 
Tarter et al. 2007). The excess absorption around the 
planet could produce a similar transit light-curve sig¬ 
nature; however mirrors are equally opaque at all visi¬ 
ble wavelengths, while atmospheric transmission exhibits 
wavelength-dependent features. Extensive optical/UV 
transit spectroscopy observations of atmospheric escape 
from hot jupiters HD 209458b reveal that escaping ma¬ 
terial contai ns elements such as H, C, O, and Na (see 
for example Ehrenreichl 2011[ |Linsky et J 2010 Char- 
bonneau et al.||2002 ), effects not expected for engineered 


structures. Lyman-a transit observations of HD 189733b 
show temporal variability in the evap orating atmosphere 
(Lecavelier des Etangs et al. 2012), whereas a stable 
mirror Meet will produce no such fluctu ations. As a 
related example, 
sit depth variabi. 


Rappaport et al. (2012) explain tran- 
ity tor KIC 12557548 b with a trail¬ 


ing cloud of material escaping a disintegrating super¬ 
mercury; this phase-folded transit light curve exhibits 
substantial asymmetries. Here transit depth should de¬ 
pend on wavelength, altho ugh multi-wavele ngth studies 
have not yet detected this dCroll et al.|2014|). In detailed 
modeling by [Brogi et al. PT2| a pre-transit flux ex¬ 
cess is due to forward scattering by dust, another feature 
not anticipated with a mirror fleet designed for optimal 
e fficiency. 


Kipping (2014) finds that bipolar residual signatures 
(somewhat different in functional form from those dis¬ 
cussed in Section 0 can result from uncorrected low- 
amplitude transit timing variations or from uncorrected 
variations in transit duration. In principle, through trial 
and error it should be possible to offset variations in ei¬ 
ther transit timing or transit duration and sharpen the 
transit entry and exit for the folded light curve, mak¬ 
ing residuals consistent with those for an isolated planet. 
However if the bipolar residuals are due to a mirror 
fleet, they will not be improved by any small adjust¬ 
ments to the timing or duration of each eclipse. Each 
subsesquent transit would make the difference more pro¬ 
nounced. Substantial technological improvements should 
eventually yield sensitivity sufficient for detailed analy¬ 
sis of individual transit light curves, in which case the 
mirror fleet would be clearly distinguishable from these 
effects. 


6. CONCLUSIONS 

An orbiting fleet of mirrors poses significant but pos¬ 
sibly not insurmountable environmental, economic, and 
engineering challenges, so could plausibly be utilized to 
redirect starlight onto a planet, perhaps to illuminate 
its dark side. Differences in the time-span of planetary 
transits compared with planetary eclipses could be one 
indication of the presence of such structures, which would 
also affect the transit light curve. 

We modeled transit light curves resulting from such 
fleets; compared with an isolated-planet transit of the 
same depth, the transit with orbiting mirrors begins ear¬ 
lier and reaches its minimum more gradually, with dif¬ 
ferences most apparent near transit entry and exit. The 
characteristic signature is not plausibly reproduced by 
alterations in limb-darkening due to abundance varia¬ 
tions. Regardless of stellar type, the difference is typi¬ 
cally 4-6% of the transit depth for mirrors extending 3 
planetary radii beyond a planet with Rp = 2i?Earth- Al¬ 
terations to planetary scattered light and thermal emis¬ 
sion caused by the mirrors are negligible compared to 
the mirrors’ effects on the light curve. For smaller stars, 
mirror fleets alter transit depth by a larger fraction of 
the total starlight. The difference is of order 10“"^/star 
for an MO star, roughly a factor of 2 smaller for a G2 
star, and several times larger for an M5 star. Effects on 
light curves are most apparent for large planets and mir¬ 
ror fleets of large radial extent. For orbital inclinations 
|j| less than 90° (zmin ^ 0), mirror effects are evident 
during a greater fraction of the transit time-scale, and 
produce inconsistencies between transit depth and time. 
Variations in mirror fleet distribution produce very small 
alterations in the light curve, so a constant-absorptance 
model sufficiently captures the main features of such an 
engineering structure. 

The mirror installations considered here are unde- 
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tectable in the existing Kepler data. However, with ap¬ 
propriate analysis, they could be revealed in JWST ob¬ 
servations, were they present. Natural systems, such as 
rings or an evaporating atmosphere, might produce light 
curves with a similar signature, but the effects will vary 
with wavelength quite differently than the mirror ffeet 
case. With sufficient follow-up observations and sensi¬ 
tivity, a mirror fleet should be distinguishable from these 
effects. 

These models likely do not accurately represent a real¬ 
istic fleet of orbiting mirrors each adjusting its attitude to 
maintain a desired illumination pattern. We anticipate 
this effect is of similar scale to differences in mirror fleet 
geometry. Future work will include detailed modeling 
of the absorptance / transmittance of plausible mirror 
fleets, incorporating consideration of the optimal range of 
mirror orbit size and eccentricity to provide appropriate 
illumination, and variations of mirror angle throughout 
the orbit. A similar effort needs to be made on whether 
stable mirror orbits are achievable considering tidal ef¬ 
fects of the star and radiative effects of the illumination 
on the thin mirrors, as well as the feasibility of propulsion 
systems capable of sustaining orbits under these condi¬ 
tions. 


This research has made use of the Exoplanet Or¬ 
bit Database and Exoplanet Data Explorer at exoplan¬ 
ets.org and NASA’s Astrophysics Data System. Dur¬ 
ing this work EJK was supported in part by NASA 
grant NNX09AN69G and donations from the Friends of 
SETI@home (setiathome.berkeley.edu). EJK would like 
to thank SMS for allowing him the honor of first au¬ 
thorship for conceiving the original idea and describing 
the basic analysis technique, despite the fact that SMS 
executed most of the work and did most of the writing. 
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